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Abstract. Entamoeba histolytica and — Spironucleus 
barkhanus have senes that encode short iron-dependent 
hydrogenases (Fe-hydrogenases), even though these protists 
lack hydrogenosomes. To understand better the biochemis- 
try of the protist Fe-hydrogenases, we prepared a recombi- 
nant E. histolytica short Fe-hydrogenase and measured its 
activity in vitro. A Giardia lamblia gene encoding a short 
Fe-hydrogenase was identified from shotgun genomic se- 
quences, and RT-PCR showed that cultured entamoebas and 
giardias transcribe short Fe-hydrogenase mRNAs. A second 
E. histolytica gene. which encoded a long Fe-hydrogenase, 
was identified [rom shotgun genomic sequences. Phyloge- 
netic analyses suggested that the short Fe-hydrogenase 
genes of entamoeba and diplomonads share a common 
ancestor, while the long Fe-hydrogenase gene of entamoeba 
appears to have been laterally transferred from a bacterium. 
These results are discussed in the context of competing 
ideas for the origins of genes encoding fermentation en- 
zymes of these protists. 


Introduction 


One of the great recent discoveries in the cell biology of 


protists 1s that the hydrogenosome of Trichomonas vagina- 
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lis, cause of vaginitis, is a modified mitochondrion, in which 
the enzymes of oxidative phosphorylation have been re- 
placed by fermentation enzymes that produce hydrogen gas 
(Müller. 1993, 1998: Bui et al., 1996; Horner et al., 1996; 
Andersson and Kurland. 1999; Rotte et al., 2000). Proof of 
this idea includes the presence of mitochondrion-like chap- 
erones and a mitochondrion-lke ATP/ADP transporter 
within the hydrogenosome, as well as organelle-targeüng 
sequences at the N-termini of hydrogenosomal proteins that 
are encoded in the nucleus (Johnson et al., 1990: Hrdy and 
Müller. 19952, b; Bui et al., 1996; Bui and Johnson, 1996; 
Horner et al., 1996: Bradley et al., 1997; Dyall et al., 2000). 

The common origin of hydrogenosomes and mitochon- 
dria is included in a new biochemical explanation For the 
origin of mitochondria, called the hydrogen hypothesis 
(Martin and Müller, 1998). The hydrogen hypothesis, à 
revision of the widely accepted endosymbiont hypothesis 
(Gray et al., 1999), suggests that the a-proteobacterium, 
which became the mitochondrion, was a facultative anaer- 
obe that was selected for its ability to produce hydrogen in 
a methanogenic archaeal host. Consistent with this idea. 
multiple hydrogenosomal fermentation enzymes—includ- 
ing ferredoxin, succinyl-CoA synthetase, and malic en- 
zyme—resemble their counterparts in mitochondria (John- 
son et al., 1990; Hrdy and Müller, 1995b). One alternative 
hypothesis suggests that the mitochondrial endosymbiont 
was selected for its ability to consume oxygen and thus 
protect the proto-eukaryote from oxidative damage 
(Andersson and Kurland, 1999). Another alternative hy- 
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pothesis suggests that the endosymbiosis was based upon 
cycling of sulfur (Searcy, 1992). 

An iron-dependent hydrogenase (Fe-hydrogenase), which 
gives the hydrogenosome its name, transfers electrons from 
reduced ferredoxin to two protons to make hydrogen gas 
(Lindmark and Müller, 1973; Müller, 1993: Payne et adi- 
1993; Bui and Johnson, 1996; Horner er al, 2000). Fe- 
hydrogenase activity. which is detected using hydrogen gas 
and methyl viologen us a reporter, has been shown in 
extracts of eukaryotes with hydrogenosomes or plastids 
(Payne et al., 1993; Wunschiers et al.. 2001). Fe-hydroge- 
nases, which are also present in strictly anaerobic bacteria, 
are designated as long or short depending upon the number 
of N-half ferredoxin-like domains that are adjacent to a 
conserved C-half hydrogenase domain (Cammark, 1992). 
The structures of a short Fe-hydrogenase of Desulfovibrio 
desulfuricans and a long Fe-hydrogenase of Clostridium 
pasteurianum have been solved, und each revealed two 
ferredoxin-like [4Fe-4S] iron-sulfur centers and a hydroge- 
nase active-site composed of a [4Fe-4S] center bridged to a 
[2Fe] eluster (Peters et al.. 1998: Nicolet et al., 1999). 

The reduced ferredoxin is produced within the hydro- 
genosome hy pyruvate:ferredoxin oxidoreductase (PFOR), 
which decarboxylates pyruvate to acetyl-CoA and CO, 
(Lindmark and Müller, 1973; Müller, 1993; Payne et al., 
1993; Hrdy and Müller, 1995a; Horner et al., 1999). It has 
been difficult to determine the phylogeny of the genes that 
encode Fe-hydrogenase and PFOR. because these iron-sul- 
fur proteins are absent from mitochondria, a-proteobacteria, 
and most other eubacteria. The pfor genes of T. vaginalis, 
Entamoeba histolytica (cause of amoebic dysentery), and 
Giardia lamblia (a cause of diarrhea) appear to share a 
common ancestry, although no bacterial donor has been 
identified, and a gene encoding a second keto-acid oxi- 
doreductase of G. lamblia appears to have a distinct ances- 
try (Reeves, 1984; Hrdy and Müller. 1995a; Rosenthal et 
aero Biown er al 1995: Müller. 1995: MOmenD cpa. 


— 


1999; Huston and Petri, 2001). The fe-hivdrogenase genes of 


T. vaginalis, E. histolytica, and Spironucleus barkhanus (a 
diplomonad similar to G. lamblia). each of which encodes a 
short Fe-hydrogenase, also appear to share a common an- 
cestry, although no bacterial donor was identified (Horner et 
al., 2000; van Hoek et al.. 2000). In contrast. the long 
fe-hydrogenase gene of the microaerophilic ciliate Nycto- 
therus ovalis appears to have a distinct ancestry. 

The Fe-hydrogenase results are surprising for three rea- 
sons (Horner et al., 2000). First. the fermentation enzymes 
of E. histolytica and G. lamblia ave present in the cytosol 
rather than in an organelle (Reeves, 1984: Müller, 1993; 
Brown et al., 1998; Mai et al., 1999; Ghosh er al., 2000). 
Entamoebas have a mitochondnon-derived organelle called 
the erypton or mitosome, while the giardial gene encoding 
a 60-kDa chaperonin appears to be endosymbiont-derived, 
but no organelle has been identified (Clark and Roger, 1995; 


Roger et al, 1995; Mai et ala 1ID99:"Tovareer ai, 1999: 
Ghosh er al., 2000). Second, entamoebas and diplomonads 
have long been thought not to produce hydrogen gas in 
culture, although a recent report suggests giardias may 
produce hydrogen under anaerobic conditions (Lindmark 
and Müller. 1973; Reeves. 1984; Brown et al., 1998; Mül- 
ler, 1998: Lloyd and Harris, 2002). Third. short fe-Aydro- 
genase genes of E. histolytica and S. barkhanus appear to 
share a common ancestor, even though these two protists are 
not closely related to each other in phylogenies drawn with 
rRNA or protein sequences (Sogin and Silberman, 1998; 
Horner et al., 2000: Bapteste et af, 2002). The latter result 
suggested the possibility of lateral transfer of fe-/ivdroge- 
nase genes between these protists. Recent. phylogenetic 
studies suggest that numerous genes encoding fermentation 
enzymes and other proteins of E. histolytica, G. lamblia, 
and 7. vaginalis may have been laterally transferred from 
prokaryotes (Rosenthal et al., 1997; Doolittle, 1998, 1999; 
de Koning et al., 2000: Field et al., 2000: Nixon et al., 
2002). 

With the goal of understanding better the biochemistry 
and evolution of Fe-hydrogenases of entamoebas and 
diplomonads, we performed the following studies. (1) The 
E. histolytica short Fe-hydrogenase | was expressed as a 
glutathione-S-transferase (GST) fusion-protein in E. coli, 
and its activity was measured 74 vitro. (2) The G. lamblia 


fe-hydrogenase gene was identified from shotgun genomic 


sequences, and mRNAs encoding short Fe-hydrogenases 
were detected in cultured giardias and entamoebas. (3) An 
E. histolytica gene encoding a long Fe-hydrogenase 2 was 
identified from genomic sequences and compared with other 
long Fe-hydrogenases. (4) Phylogenetic analyses were re- 
peated with the addition of Fe-hydrogenases of G. lamblia, 
E. histolytica, the green alga Chlamydomonas reinhardtii, 
and the eubacterium Alegasphaera elsdenii. 


Materials and Methods 


Cloning of the G. lamblia fe-hydrogenase and E. histo- 
lytica genes and identification of mRNAs encoding Fe- 
hydrogenase from cultured entamoebas and giardias. An E. 
histolytica EST (GenBank ABOO2772). which encodes the 
N-terminus of a putative short Fe-hydrogenase 1, was iden- 
tified from GenBank using BLASTP and an amebic 2[+Fe- 
4S] ferredoxin sequence (Altschul et al., 1997; Nixon et al, 
2002). The 3’ end of the E. histolytica fe-hydrogenase 1 
gene was isolated using 3’ RACE (FirstChoice RLM-RACE 
kit, Ambion lnc., Austin, Texas), and the entire entamoebic 


fe-hydrogenase 1. gene was cloned, sequenced on both 


strands, and deposited in GenBank under accession number 
AAGO09783. A second Entamoeba fe-lvdrogenase gene, 
which encodes a long Fe-hydrogenase, was identified from 
assemblies of E. histolytica genome sequences at The Insti- 
tute for Genomic Research. The E. histolytica fe-hyvdroge- 
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nase 2 gene was deposited in GenBank under accession 
number AY 172963. 

A shotgun clone from the G. lamblia genome sequencing 
project, which contained the 5' end of a putative fe-livdro- 
genase gene, was used to identify the entire fe-/vdrogenase 
gene from a G. lamblia genomic DNA library made in 
Lambda Zap (McArthur et al., 2000). The G. lamblia fe- 
hydrogenase gene was sequenced on both strands and de- 
posited in GenBank under accession number AAK28337. 
The N-termini of predicted entamoebic and giardial Fe- 
hydrogenases were examined with MITOP to determine 
whether organelle-targeting sequences might be present 
(Claros and Vincens, 1996). 

Total RNA was prepared by lysing cultured entamoebas 
and giardias in a guanidinium isothiocyanate solution and 
by centrifuging the lysate through a cesium chloride gradi- 
ent (Choczynski and Sacchi, 1987). Reverse-transcriptase 
and polymerase chain reaction (RT-PCR) were performed 
with these RNAs and with primers specific for the E. 
histolytica and G. lamblia genes encoding short Fe-hydro- 
genase, malic enzyme, alcohol dehydrogenase E (ADHE), 
and ferredoxin | (entamoebas only) (Rosenthal ef al.. 1997; 
Nixon ef al, 2002). For negative controls, PCR was per- 
formed without RT, and RT-PCR products were identified 
on agarose gels. 

Expression of a recombinant short E. histolytica Fe- 
livdrogenase and measurement of its activity. A recombi- 
nant glutathione-S-transferase (GST) fusion-protein. con- 
taining a short Fe-hydrogenase at its C-terminus was made 
by cloning the £. histolytica Fe-hydrogenase | coding re- 
gion into the pGEX-6T vector (Smith and Johnson, 1988). 
The GST-Fe-hydrogenase | construct was transfected into 
Escherichia coli strain BL21, whieh was grown anaerobi- 
cally and induced with isopropyl B-D-thiogalactopyrano- 
side (IPTG). Bacteria were lysed by freezing and thawing. 
and the hydrogenase activities of the supernatant and of 
purified GST-Fe-hydrogenase | were determined. As a neg- 
ative control, the activity of bacterial lysate expressing a 
GST-chitinase fusion protein was measured. In a septum- 
sealed cuvette, Fe-hydrogenase ] activity was exanuned by 
incubation with 10 mM methyl viologen in 50 mM Tris-HCl 
(pH 8.0) and 5 mM dithioreitol, which had been bubbled 
with hydrogen for 15 min (Payne et al., 1993). Reduction of 
methyl viologen was monitored at 600 nm. A negative 
control was bubbled with nitrogen instead of hydrogen. 

We also measured Fe-hydrogenase activity from lysates 
of cultures of non-transformed 7. vaginalis and trom E. 
histolvtica that had been transformed with à plasmid con- 
taining the entamoebic fe-/ivdrogenase / gene under an 
actin gene promoter (Ghosh er al., 2000). This construct 
was previously used to demonstrate the cytosolic location of 
the entamoebic Fe-hydrogenase 1. 

Phylogenetic analyses. Amino acid sequences homolo- 
gous to the G. lamblia and E. histolytica Fe-hydrogenases, 


which were identified in the nonredundant and unfinished 
microbial databases of GenBank using BLASTP and 
TBLASTN, respectively, were aligned using CLUSTALW 
(Thompson et al., 1994; Altschul er al., 1997). The align- 
ment was adjusted manually using the SEQL AB program 
(Genetics Computer Group, Madison. Wisconsin). Regions 
that could not be unambiguously aligned. were excluded, 
leaving 351 aligned amino acid positions for use in phylo- 
genetic analyses. Phylogenetic analyses were performed 
using distance and parsimony methods by the computer 
programs TREE-PUZZLE (Strimmer and von Haeseler, 
1996) and PHYLIP (Felsenstein, 1989). Pairwise distances 
were computed using TREE-PUZZLE under the Dayhoff 
model (Dayhoff et al., 1978), with the inclusion of observed 
amino acid frequencies, estimated proportion of invariant 
sites, and estimation of among-site variation for the remain- 
ing sites according to a gamma distribution (four discrete 
categories). The optimal tree was inferred using the Fitch- 
Margohash algorithm (Fitch and Margoliash. 1967), using 
the Fitch program, with global rearrangements and 100 
random-addition replicates. Bootstrap values were obtained, 
using the 100 resampled datasets. under the same model 
using the PUZZLEBOOT program (available at http://www. 
tree-puzzle.de). The optimal tree under parsimony and re- 
lated bootstrap values were determined using PHYLIP's 
PROTPARS program for parsimony. 


Results and Discussion 


A recombinant short Fe-hydrogenase | of E. histolytica is 
active. The amino-terminus of the E. histolytica Fe-hydro- 
genase | did not contain an organelle-targeting sequence 
(Claros and Vincens, 1996: Horner et al., 2000). Indeed, an 
epitope-tagged Fe-hydrogenase | is present within the cy- 
tosol of transfected E. histolytica (Ghosh et al., 2000). The 
activity of the entamoebic Fe-hydrogenase ] was measured 
in lysates of E. coli, which were expressing a GST-ent- 
amoebic Fe-hydrogenase | fusion-enzyme (Table 1) (Smith 
and Johnson, 1988). The activity of the recombinant ent- 
amoebic Fe-hydrogenase ] was present when hydrogen was 
bubbled into the medium but absent when nitrogen (nega- 
tive control) was bubbled into the medium (Payne er al., 
1993). There was no Fe-hydrogenase activity in control E. 
coli, which were overexpressing a GST-chitinase fusion- 
protein. Recombinant entamoebic Fe-hydrogenase 1, which 
was purified on glutathione-agarose beads, had a decreased 
specific activity (data not shown). This was likely caused by 
exposure to oxygen during the purification procedure that 
probably inactivated the Fe-hydrogenase | iron-sulfur cen- 
ters (CammsebeaPsvne ef al. 1993: Horner ef al. 
2000). The specific activity of the GST-Fe-hydrogenase | 
fusion enzyme was about 10 times that of the native ent- 
amoebic Fe-hydrogenase ] overexpressed in transfected E. 
histolvtica (Table 1). Interestingly, the specific activity of 
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Table t 
Activities of Entamoeba histolytica and Trichomonas vaginalis 


Fe-livdrogenases 


Hydrogenase activily 


Sample (nmol/min/mg of protein)" 


Bacteria transformed with E. histolytica 


Fe-hydrogenase 36 (2) 
Transfected E. histolytica with Fe- 

hydrogenase 300 (3) 
Trichomonas vaginalis 167 + 32 (9) 
T. vaginalis + 0.47 mg E. histolytica 

lysate” 114 (2) 


" Averages +/— standard deviations, where possible, Number of deter- 
minations in parentheses. 

^ Calculated A, = 0.56 mg (amount of E. histolytica lysate m mg of 
protein to cause 50% reduction of T. vaginalis hydrogenase activity). 


the T. vaginalis Fe-hydrogenase was greater than that of the 
recombinant entamoebic GST-entamoebic Fe-hydrogenase 
and was inhibited by a lysate of non-transfected entamoebas 
(Table 1). This may explain why it was difficult to detect 
Fe-hydrogenase activity in lysates of nontransfected ent- 
amoebas, even though Fe-hydrogenase ] mRNAs were 
identified from them by RT-PCR (next section). 

Cultured entamoebas and giardias express n RNAs en- 
coding short Fe-hydrogenases. We isolated an fe-hydroge- 
nase gene of G. lamblia, because we have frequently com- 
pared the fermentation enzymes of this diplomonad with 
Ihose of E. histolytica (Rosenthal et al.. 1997; Field er al.. 
2000; Nixon et al., 2002). A search of the contigs predicted 
from the G. lumblia shogun sequences suggested that this 
gene, which predicts a short Fe-hydrogenase, is the only 
hydrogenase gene present within the giardial genome. Like 
the entamoebic. Fe-hydrogenase 1, the predicted. giardial 
Fe-hydrogenase lacked an N-terminal organelle-targeting 
sequence and had two ferredoxin-like iron-sulfur centers 
and a hydrogenase iron-sulfur center like those present in 
the short Fe-hydrogenases of T. vaginalis, Desulfovibrio sp., 
and Clostridia sp. (Cammark, 1992; Thompson et al., 1994; 
Bui and Johnson, 1996; Horner et al.. 1996; Nicolet et al., 
1999), RT-PCR showed that cultured entamoebas and giar- 
dias contain mRNAs, which encode short Fe-hydrogenases 
(Fig. 1A. B). Negative controls without RT showed that the 
RT-PCR was not amplifying DNA from the extracts of 
cultured entamoebas and giardias. Because the giardial con- 
tigs predicted only one Fe-hydrogenase, which is expressed, 
it is likely that the hydrogenase activity recently detected in 
cultures of giardias derives Irom this enzyme (Lloyd and 
Harris. 2002). In contrast, entamoebas appear to have a 
second long hydrogenase (see next section), so if entamoe- 
bic hydrogenase activity is present, it might derive from one 
or more enzymes. These results suggest the possibility that 
entamoebas and giardias use protons as electron acceptors 
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Figure t. Agarose gels of ethidium-stained RT-PCR products from 
entamoehic and giardial mRNAs. Images are reversed for clarity of 
reproduction. (A) RT-PCR of amochic mRNAs encoding malic enzyme 
(lane 1). Fe-hydrogenase (lane 2). alcohol dehydrogenase E (lane 3). 
and ferredoxin (lane 4). Size markers are shown in lane 5. Lanes 6-9 
are the negative controls for malic enzyme. hydrogenase. ADHE, 
and ferredoxin, respectively, (B) RT-PCR of giardial mRNAs encod- 
ing ADHE (lane 3), Fe-hydrogenase (lane 4). and mahe enzyme (lane 
5). A negative control (no RT) for Fe-hydrogenase is shown in lane 6. 
A positive control for Fe-hydrogenase. using Giardia lamblia WB 
strain DNA, is shown in lane 2. Size markers are shown in lane t. 
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MSTOLTPLRNKIISEVVKCFKSGRFIEDIDKLPTILTDGDGWKPTSKFVHSREQEEGIYR 
DLKREDBVRIARDOFEGKLOEGVHYTPREMVPRN.STEI.RCCIBEHDR,.EIMR 
VRHKLLAKLVNLWKENKLTNEIDRLPIELSPRR.SRPLGRCCIHKER..AVYK 


EKVLSVLGF.VDGEYDDITPLHVYAOKALERT.SLHEPVFGISQKGCNKCHFNGYFVTQA 
HRVIARLGCSLENYDEEKT.LAOFAKEALERE.KPTWPMLTVLDEACNACVKSKYMITNA 
YKLFPLLGFDMTDETDELTSLSEYAROALERKNKOKENILCVIDEACSSCVOVNYEVTNL 


* * 


CEGCTSRPCSVNCPRKCISFGEDGRAVINQNNCIKCGRCYKFCPYGAIISKSVPCVKACP 
COACVARPCMMNCPKTAIAIS.GGRARIDEEKCINCGICLKNCPYHAVIKIPVPCEEACP 
CRGCVARSCYMNCPKDAIRFRKNGOAKIDHDACISCGKCHQSCPYHAIVFIPVPCEEACP 
* * * * x x x x * * 


CGAMLDSPEGVKTIDFEKCINCGGCMRACPFGAILPRSNLIDVLK.ILPTKKVVACPAPS 

VGAISKDENGKERIDYHRKCI FCGNCMRECP FGAMMDKGOIVDVI KHLMSGKKVSALYAPA 

VKAISKDENGIEHIDESKCI YCGKCLNACPFGAT FEISQAFDVLEGIRSGEKMIAIPAPS 
X X PA X 


IAAHFGKYDLALVSGGLIOVGFTSVEDVSYGADLCALNEAKEFEERIVKNKKDFMTTSCC 
VAAOF.KAVPGOLESALKKAGFNKVWEVAIGADITADREASEFEERMEHGHI.LMTTSCC 
ILGOF.NTSIEAVYGALROMGFADVVEVAOGAMDTVSHEAAELKEKLEEGQP.FMTTSCC 

O 


PAYINAINKHMPELKENVSHTPTPMHFATOAVKDRDQETVTVFIGPCNAKRWETLODSTT 
PAYVRAVKKHVPALVPCISDTRSPMHYTAELAKKEDPDCVTVFIGPCLAKRREGLEDEFV 
PSYIELVNKHIPGMKPYVSSTGSPMYYAARIAKERHPDAKIVFIGPCVAKRKEARRDECV 


O 


DYCLTFDEIFGLFEGSGIDLSKVOPYTFVDKAHKEGKIFAVSGGVASAVASLLPKEVPDG 
DYVLSIEELGALLTAKEIDISKEEALPGKITPTSSGRGFAASGGVAEAVRVRL.KKPEN. 
DYILTFEEMASIFEGLDIOLEOTOPFSVLYTSVREAHGEAQAGGVMGAIKAYLGEEAKK. 


ENEPTITDDGFSOENERRERNERENI. vs TGNLVEVMVCEGGCAYGPGCPGLNTP 
. LRPVLINGLNKEGMKOLAS YGKIOSGELPHDSSTPNLVEVMSCEGGCIGGP 
SESATOVSDLNKRENIGLLRAAAKTG......-. KAQGOFIEVMACEGGCISGP 


O O 
ATSAKIKIAVDKMEAHPEGRWVGLPNSQIKPIKVEN 504 


Alignment of the predicted Entamoeba histolytica long Fe-hydrogenase 2 (Eh2) with predicted 


long Fe-hydrogenases of Treponema denticola (Td) und Bacteroides fragilis (Bf). Conserved Cys residues, 
which are shaded, include those that coordinate putative [4Fe-4S] iron-sulfur centers (marked with x°s) and Ihose 
that coordinate putative hydrogenase iron-sulfur centers (marked with o's). Other conserved Cys residues, which 
may be involved in coordinating iron-sulfur centers. are marked with asterisks. Amino acids at the beginning and 


end of the conserved Fe-hydrogenase domain are underlined. 


when the organisms are growing under strictly anaerobic 
conditions in the bowel lumen (Brown et al., 1998; Huston 
and Petri, 2001; Lloyd and Harris, 2002). 

E. histolytica has a hydrogenase 2 gene encoding a long 
Fe-hydrogenase. The assembhes of the shotgun sequences 
of the E. histolytica genome predicted a long Fe-hydroge- 
nase 2 (Fig. 2) in addition to the short Fe-hydrogenase 1. 
The entamoebic Fe-hydrogenase 2 was 504 amino acids 
long and had an N-terminal sequenee, which included pos- 
itively charged Lys and Arg that are often present at the 
N-termini of organellar proteins (Claros and Vincens, 
1996). In addition, the N-terminus of Fe-hydrogenase 2 


— 


contained Ser and Leu residues, which are present at the 
N-termini of crypton and hydrogenosomal proteins (Bui ef 
al.. 1996; Materai 1999) However, in the absence of 
experimental evidence, we cannot be sure that the entamoe- 
bic long Fe-hydrogenase is targeted to the erypton. 

The entamoebic Fe-hydrogenase 2 was much more sim- 
ilar (3856 amino acid identities) to predicted long Fe- 
hydrogenases of Bacteroides fragilis and Treponema den- 
ticola than to short Fe-hydrogenases of entamoebas, giar- 
dias, trichomonads, and other anaerobic bacteria (<28% 
amino acid identities; Fig. 2). The entamoebic Fe-hydroge- 
nase 2 and the predicted long Fe-hydrogenases of B. fragilis 
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and T. denticola each contained Cys residues that likely 
coordinate two ferredoxin-like [4Fe-4S] iron-sulfur centers 
(marked with x's in Fig. 2) and hydrogenase iron-sulfur 
centers (marked with o's), which have previously been 
identified in structures of short and long Fe-hydrogenases 
(Peters et al. 1998; Nicolet et a/., 1999). In addition, the 
predicted entamoebic Fe-hydrogenase 2 had eight other 
N-terminal Cys residues, which aligned with those of the 
bacteroides and treponema long Fe-hydrogenases (marked 
with asterisks). Although these Cys residues probably co- 
ordinate other iron-sulfur centers, they remain unidentified, 
because they do not align with the N-terminal iron-sulfur 
centers of the long Fe-hydrogenase of C. pasteurtanium, 
which has been crystallized (Peters eral., 1998), 

The entamoebic and giardial short fe-hydrogenase | 
genes appear to share a common ancestry, while the ent- 
amoebic long te-hydrogenase 2 gene appears to have been 
laterally transferred from a prokaryote. Phylogenetic trees 
of Fe-hydrogenases from eubacteria and eukaryotes are 
star-shaped and contain few basal nodes that are strongly 
supported (Fig. 3). This result suggests that the Fe-hydro- 
genases are widely divergent and that little phylogenetic 
signal remains. For example. Fe-hydrogenases of closely 
related eukaryotes—either trichomonads, green algae 
(Chlamydomonas reinhardtii, Scenedesmus obliquus, and 
Chlorella fusca), or chytrid fungi (Prromyces sp. and Neo- 
callimastix frontalis)— each grouped together, but Fe-hydrog- 
enases of unrelated eukaryotes did not group together. 1n 
particular, our analysis does not support recent conclusions 
that hydrogenases of trichomonads are monophyletic with 
those of chytrid fungi (Voncken et al., 2002) or with those 
of E. histolytica and S. barkhanus (Horner et al., 2000). 

The short fe-/aivdrogenase genes of G. lamblia, S. barkla- 
nus, and E. histolytica appear to share a most recent com- 
mon ancestry, although a particular bacterial donor was not 
identified. Remarkably, the short Fe-hydrogenase of G. lam- 
blia was more similar to that of E. histolytica than to that of 


S. barkhanus. Because G. lamblia and S. barkhanus are 
diplomonads, which share à recent common ancestor in 
phylogenetic trees of rRNA and proteins (Sogin and Silber- 
man, 1998), a possible explanation of these results is that 
the E. histolytica fe-hydrogenase gene was laterally trans- 
ferred from a diplomonad (Rosenthal et al., 1997; Doolittle, 
1998, 1999; Müller, 1998; de Koning et al., 2000; Field et 
al., 2000: Nixon et al., 2002). This lateral gene transfer 
would not have occurred recently, because the Fe-hydroge- 
nases of entamoebas and giardias showed only a 40% amino 
acid identity with each other, and each fe-/ivdrogenase gene 
has the codon usage of its host. Alternatively, the 
diplomonad-F. histolytica sub-clade could be incorrectly 
rooted by the long branch connecting it to the remainder of 
the tree. 

The common ancestry of genes encoding the E. histo- 
lvtica long Fe-hydrogenase 2 and those of B. fragilis and T. 
denticola is strongly supported. This appears then to be an 
example of lateral gene transfer, as Entamoeba is not a close 
relative of either of these eubacteria (Rosenthal et al., 1997; 
Doolittle, 1998, 1999: Müller, 1998; de Koning et al., 2000; 
Field et al., 2000; Nixon et al.. 2002). There was weak 
support for the pairing of Pe-hydrogenases of the ciliate N. 
ovalis and Desulfovibrio sp., as has been previously noted 
(Horner et al., 2000; Voncken et al., 2002). This suggests 
that the ciliate hydrogenase was derived by lateral gene 
transfer, but does not prove it. 


Conclusions 


This is the first time that an Fe-hydrogenase from a protist 
has been expressed as a GST-fusion protein in bacteria. This 
is also the first time that an fe-/iydrogenase gene (encoding 
the long hydrogenase of entamoebas) has been inferred to 
have been laterally transferred from a bacterium, although 
numerous genes encoding fermentation enzymes (e.g., al- 
cohol dehydrogenases, malic enzyme, and acetyl-CoA syn- 


Figure 3. Phylogenetic relationships of Fe-hydrogenases. inferred using a distance matrix generated by the 
Dayhoff+1+P model and the Fitch-Margoliash algorithm. Bootstrap values obtained using PUZZLEBOOT and 
PROTPARS, respectively. are shown at the relevant nodes. Bootstrap values below 50% are marked with an 
asterisk if the other bootstrap value is >50%. 1f both hootstrap values are below 50%, neither is marked. The 
scale bar indicates estimated sequence divergence per unit branch length. Sequences from eukaryotes, which are 
boxed. were from Chlamydomonas reinhardtii | and 2 [accession # 16945126 and 18026272]: Chlorella fusca 
[21732235]; Entamoeba histolytica | and 2 (9963974 and AY 172963]; Giardia lamblia [13506793]: Neocal- 
limastix frontalis 119547863]: Nycrotherus ovalis [4034791]; Prromyces sp. [19548180]: Scenedesmus obliquus 
] and 2 112581498 and 13311187]; Spironucleus barkhanus [11127703]; Trichomonas gallinae [19548182]; 
Trichomonas vaginalis |, 2, 3, and 4 | 19547859, 11127701, 1171117, and 1345094]. Eubacterial sequences were 
from Bacteroides fragilis [unfinished microbial database]: Clostridium acetobutylicum [15896476]; Clostridium 
pasteurianum [557064]; Clostridium perfringens 1,2. and 3 (18311557. 18309258, and 18311328]; Clostridium 
saccharobutylicum [458597]: Clostridium tlhermocellum [3927278]. Desulfovibrio fructosovorans [1914864]: 
Desulfovibrio desulfuricans | and 2 [4930044 and 13022069]: Desulfovibrio vulgaris | and 2 [97381 and 66319, 
respectively], Enbacterium acidominophilum |14250935), Megasphacra elsdenii [6650985]; Ruminococcus 
albus [unfinished microbial database]: Thermoanaerobacter tengcongensis | and 2 [20807184 and 20515894]: 
Thermotoga maritima 1, 2, and 3 (15644177, 7433127. and 4981985]; and Treponema denticola [unfinished 


microbial database]. 
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thases) appear io have been laterally transferred from pro- 
karyotes to amoebas and gtardias (Rosenthal er al.. 1997: 
Field et al., 2000; Nixon et al.. 2002). Although the evi- 
dence is weak, this may also be the first time that a gene 
(encoding the short hydrogenase of entamoebas) has been 
inferred to have been laterally transferred. from another 
prolist. Because the hypothesized lateral gene iransfer 
would probably have occurred atier the acquisition of the 
fe-hydrogenase gene by the diplomonad lineage, this par- 
ticular result does not disprove ihe hydrogen hypothesis 
(Martin and Müller, 1998). However, the failure to demon- 
sirate that the eukaryotie Fe-hydrogenases share à common 
ancestry, or to identify an a-proteobacterial donor for these 
eukaryotic fe-hydrogenase. genes (Horner et al.. 2000). 
dampens our enthusiasm for the hydrogen hypothesis. 
These results suggest that the mitochondrial endosymbiont 
was selecied for a properiy other than hydrogen production 
(¢.g., its ability to consume oxygen) (Andersson and Kur- 
land, 1999) and that ihe presence of Fe-hydrogenases and 
other fermentation enzymes of microaerophilic eukaryotes 
may reflect a secondary adaptation to their anaerobic envi- 
ronment (Rosenthal et al., 1997; Doolittle, 1998, 1999; de 
Koning et al.. 2000; Field et al.. 2000; Lloyd and Harris, 
2002 Nixon eral DU? 
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